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Electron-phonon heat exchange in layered nano-systems 


D. V. Angheh and S. Cojocaru^ 

^Horia Hulubei National Institute For Physics And Nuclear Engineering, RO-077125, Mdgurele, Romania 

We analyze the heat power P between electrons and phonons in thin metallic films deposited 
on free-standing dielectric membranes in a temperature range in which the phonon gas has a 
quasi two-dimensional distribution. The quantization of the electrons wavenumbers in the direc¬ 
tion perpendicular to the film surfaces lead to the formation of quasi two-dimensional electronic 
sub-bands. The electron-phonon coupling is treated in the deformation potential model and, if 
we denote by Te the electrons temperature and by Tph the phonons temperature, we find that 
P = P^°\Te) — P^^\Te,Tph)', is the power “emitted” by the electron system to the phonons 
and P(i) is the power “absorbed” by the electrons from the phonons. 

Due to the quantization of the electronic states, P vs (d, Te) and P vs {d,Tph) show very strong 
oscillations with d, forming sharp crests almost parallel to the temperature axes. In the valleys 
between the crests, P oc Te ® — T®^®. From valley to crest, P increases by more than one order of 
magnitude and on the crests P does not have a simple power law dependence on temperature. 

The strong modulation of P with the thickness of the film may provide a way to control the 
electron-phonon heat power and the power dissipation in thin metallic hlms. Eventually the same 
mechanism may be used to detect small variations of d or surface contamination. On the other 
hand, the surface imperfections of the metallic hlms may make it difficult to observe the oscillations 
of P with d and eventually due to averaging the effects the heat how would have a more smooth 
dependence on the thickness in real experiments. 

PACS numbers: 85.85.+j, 63.20.kd, 72.15.Jf 


I. INTRODUCTION 

In a recent paper Nquyen et al.^ reported re¬ 
markable cooling properties of normal metal-insulator- 
superconductor (NIS) tunnel junctions refrigerators, by 
reaching electronic temperatures of 30 mK or below, from 
a bath temperature of 150 mK, at a cooling power of 
40 pW. Such micro-refrigerators have great potential for 
applications, since they can be mounted directly on chips 
for cooling qubits or ultra-sensitive detectors, like micro¬ 
bolometers or micro-calorimeters. 

The principle of operation of NIS micro-refrigerators 
has been explained in several publications (e.g.l2H3) and 
consists basically in cooling of a normal metal island by 
evacuating the “hot” electrons (from above the Fermi 
sea) into a superconductor while injecting “cold” elec¬ 
trons (below the Fermi sea) from another superconduc¬ 
tor using a pair of symmetrically biased NIS tunnel junc¬ 
tions. If the normal metal island is deposited on a chip, 
then it can serve as a refrigerator by cooling the chip 
through electron-phonon interaction. The efficiency of 
the electron cooling process is controlled by the bias volt¬ 
ages of the NIS junctions, whereas the success of the chip 
refrigeration strongly depends on the electron-phonon 
coupling. Moreover, due to the strong temperature de¬ 
pendence of the electric current through the junctions at 
fixed bias voltage (or the strong variation of voltage with 
temperature at fixed current) the NIS junctions can also 
serve as thermometers. Because of this, if the normal 
metal island absorbs radiation, th e device turns into a 
very sensitive radiation detector!^. 

When it works as a detector, the normal metal is¬ 
land can be kept at the nominal working temperature 


either by cooling it directly through NIS junctions (even¬ 
tually through the thermometer junctions), or indirectly 
through electron-phonon coupling to a cold substrate^Hm. 
Therefore in any situation, i.e. when the NIS junctions 
work as coolers, thermometers, or radiation detectors, 
the electron-phonon coupling plays a central role in the 
functionality of the device. 

A typical experimental setujP^ is depicted in Fig. 
and consists of a Cu film of thickness of the order of 10 nm 
deposited on a dielectric silicon-nitride (SiN^,) membrane 
of thickness of the order of 100 nm. The Cu film is the 
normal metal island and is connected to superconducting 
A1 leads through NIS tunnel junctions. When it is used 
as a radiation detector, to reach the sensitivity required 
by astronomical observations, the working temperature 
of the d evice s hould be in the range of hundreds of mK 
or belo'v\lili21iil. At such temperatures the phonon gas in 
the layered structure formed by the normal metal island 
and the supporting membrane undergoes a dimension¬ 
ality cross-over from a three-dimensional (3D) gas (at 
higher temperatures) to a quasi two-dimensional (2D) 
gas (at lower temperatures)M^l*i^. 

In the stationary regime one may assume that the elec¬ 
trons in the metallic layer have a Fermi distribution char¬ 
acterized by an effective temperature Tg, whereas the 
phonons have a Bose distribution of effective temperature 
Tph- In 3D bulk systems the heat flux between the elec¬ 
trons and phonons have been calculated by Wellstood et 
al.li^ and has been shown to vary as T® — at low tem¬ 
peratures. Such a model is not justified for our devices 
and finite size effecthave to be taken into account. 

Phenomenologically, a number of experimental stud¬ 
ies have been interpreted by assuming that the heat flux 
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FIG. 1: The schematic model of the system. On a dielectric 
membrane with parallel surfaces perpendicular on the 2 axis 
is deposited an uniform metallic him. The surfaces of the 
membrane cut the 2 axis at —1//2 and L/2 — d, wheres the 
upper surface of the him cuts the 2 axis sX L12. In the (x,y) 
plane the area of the device is A I?). The whole system is 
treated like a homogeneous elastic continuum, with the elastic 
properties and the density of silicon nitride (SiN^,). 


has a temperature dependence proportional to T® — 
where x < Qn the other hand, a theoretical in¬ 

vestigation of the surface effects for a thin metallic film 
deposited on a half-space (semi-infinite) dielectric showed 
that the value of x is actually larger than Qualita¬ 
tively, the growth of the exponent x in the lowest part of 
the measured temperature range has been later found in 
some experiments when metallic films were deposited on 
bulky substrateJi^. 

The electron scattering rate caused by interaction with 
2D (Lamb) phonon modes in a semiconductor quantum 
well (QW) has been studied in Refs. |16I21) and in a 
double heterostructure QW including the piezoelectric 
coupling in Ref. [^ . 

The electron-phonon heat transfer in monolayer and 
bilayer graphene was studied in Ref. [53] and a tempera¬ 
ture dependence of the form was found in the low 

temperature regime. For a quasi one-dimensional geom¬ 
etry (metallic nanowire) the electron-phonon power flux 
was studied theoretically in [53] and a dependence was 
obtained. It is argued that a general temperature depen¬ 
dence of the form Tg+^ — should be valid, where s 
is the smaller of dimensions of the electron and phonon 
system.l^ However, in experimental studies on quasi ID 
A1 nanowires with 65 x 90 nm^ cross-section, a better fit 
is achieved with the standard exponential x = 5.^ 

In this context it is interesting to analyze theoreti¬ 
cally the temperature dependence of the heat transfered 
between electrons and phonons in the typical experi¬ 
mental setup of Fig. [2 We assume that the metallic 
layer is made of Cu and the supporting membrane is 


silicon nitride (SiN^:). We first observ^^^ that the di¬ 
mensionality crossover of the phonon gas in a 100 nm 
thick SiNj; membrane occurs around a scaling tempera¬ 
ture Tc = cih/i^ksL) « 237 mK. Below this tempera¬ 
ture one may expect a quasi-2D behavior of the phonon 
system (rigorously, for T Tc), whereas above it (rig¬ 
orously, at T Tc) a 3D model would suffice. 

We carry out our analy sis employing a QW picture of 
the metallic filnpSiSlIllMl by taking into account the dis¬ 
cretization of the components of the electrons wavenum¬ 
bers perpendicular to the membrane’s surfaces, whereas 
the 3D Fermi gas picture of the electron system is used 
in an accompanying paper. We work in a temperature 
range T < 200 mK in which the phonon gas is quasi- 
2D and we observe that the electron-phonon heat flux 
P cannot be simply described by a single power-law de¬ 
pendence, P (xT^ — Tpi^. While in some ranges of d we 
have (X T^'^ — T^j^, in general the heat flux has a 
very strong oscillatory behavior as a function of d. At 
certain regular intervals the power flux increases sharply 
with the film thickness by at least one order of magni¬ 
tude (Fig. 1^ and in the regions of increased heat flux 
the exponent of the temperature dependence is not well 
defined (Fig. |b). 

II. ELECTRON-PHONON INTERACTION 
A. The electron gas 

The electron system is described as a gas of free 
fermions confined in the metallic film (Fig. [^. The elec¬ 
tron wavevector will be denoted by k = (kjp/cj;), where 
kj| and are the components of k perpendicular and 
parallel to z, respectively. The wavefunctions satisfy the 
Dirichlet boundary conditions on the film surfaces (at 
L/2 — d and L/2) and periodic boundary conditions in 
the {xy) plane, V'k||.n(r,t) = 
where n = 1,2,kz = im/d, and 



A is the area of the device surface. The electron energy 
is 

;j2^2 H'^k'n ^ ^ 

where mg is the electron mass, k = |k|, and fey = |k|[|. 
The electrons are distributed over energy sub-bands de¬ 
fined by n. The electron field operators are 

= X! ^k||,n(r,t)ck||,„, (3a) 

k||.n 

^^(Di) = XI ^k||.n(l':i)Ck||,n^ (3b) 

k||.n 
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where cl and Ck„ n are the electron creation and anni- 

K|| ,n *^|| 

hilation operators on the state (k||,n). 


qtQ-which are the 2 components of the wavevectors corre¬ 
sponding to the longitudinal and transversal oscillations, 
respectively-ar^^ 


1. Numerical estimations 


To better understand the physical problem, let us com¬ 
pute the magnitude of some important quantities related 
to the electron gas. The Cu Fermi energy is Ep « 7 eV; 
the highest sub-band populated at 0 K is denoted as Cnp , 
where 


np = 


^2meEp 

ttH 


d 


( 4 ) 


and [xj is the biggest integer smaller than x. For d = 
10 and 20 nm, np = A3 and 86 , respectively. The Cu 
Fermi temperature is Tp = Ep/kp ~ 81231 K, whereas 
(e 2 — ^i)lkB ~ 131 K and (£44 — e43)/kB « 3796 K for 
d = 10 nm. 


B. The phonon gas 


For simplicity, we treat the whole system-supporting 
membrane and metallic film (Fig. [^-as a single elastic 
continuum of thickness L, area A, and volume = 
LA. The elastic modes of such a systenJSEHIZIlll have 

the form where qy and ry are 

the components of the wavevector and position vector, 
respectively, parallel to the (x, y) plane. The functions 
Wq||j( 0 ) are normalized on the interval z S [—L/2, L/2], 

namely 

The elastic modes are divided into three main cate¬ 
gories: horizontal shear (h), dilatational or symmetric 
(s) with respect to the mid-plane, and flexural or anti¬ 
symmetric (a) with respect to the mid-plane. The quan¬ 
tization of the elastic modes in the z direction leads to the 
formation of phonon branches (or sub-bands), such that 
a phonon mode is identified by its symmetry a = h, s, a, 
sub-band number p = 1, 2,..., and qy; ^ represents the 
pair (a, v). 

The h modes are simple transversal modes. The free 
boundary conditions imposed at the upper and lower 
surfaces of the membrane quantize the z-component of 
the wavevector at the values qthn = np/L, n = 0,1,..., 
whereas qihn = 0. If we denote by q\\ = |qy|, the phonon 

frequency is where by Ct and ci 

we denote the transversal and longitudinal sound veloc¬ 
ities, respectively. The sound velocities are determined 
by the Lame coefficients A and fi, and the density p of 


pi A -j- 2p 

Ct = -, Cl = - 

P P 


( 5 ) 


The s and a modes are superpositions of longitudinal 
and transversal modes, oscillating in a plane perpendicu¬ 
lar to the surfaces. The quantization relations for qia and 


{ql-qLY~[tMqio.Ll2)\ ’ 

where - 1-1 and —1 exponents on the right hand side corre¬ 
spond to the a = s and a = a, respectively. The closure 
equation is secured by the Snell’s law. 


and the branch number v is determined from Eqs. (§ 
and ([7| in the limit 9 y —>■ 0. 

In the second quantization, the displacement field op¬ 
erator will be denoted by 


^,q|l 


/ h 

2pwq 5 L 


®qii€ ' 


-qii« 


X w, 




(z). 


^i(q||r||-zwq ^t) 


( 8 ) 


where Cqn^ and Oqu^ are the phonon creation and anni¬ 
hilation operators. 

We introduce the notation J = c^/c^ {J < 1/2 from 
Eq. [^. We shall write qt instead of qt^^ and qi instead 
of qp^, when this does not lead to confusions. Explicitly, 
the normalized displacement fields of the s and a waves 
ar^^ 


Wc 


= NAqt 


2 gy cos cos {qiz) + (g/ - qf) 


X cos I ^ ) cos(zgt) 


(9a) 


c^qil 


= N,q\\ 


'^qii ,a,u,x 


2qtqi cos sin (®z) -h (<?/ - gy) 

(9b) 

Naiqt 2qf^ sin sin (qiz) + (g/ - g|[) 

(9c) 


X cos ( sin(zgt) 


X sm I ^ ) sin {zqt) 


w. 


qy 


= Nr. 


n 


2qtqi sin cos (®z) - (g/ - gy) 


X sin ( ^ ) cos (zqt) 


(9d) 


where we took qy along the x axis. 

We are interested in the low temperatures limit, so gy 
is small (rigorously it should be L^y ^ 1 ) in the low¬ 
est phonon sub-band {p — 0). Under these conditions, 
qt,s,o = qt,s is real, whereas qps,o = ipi.s, ft.a.o = 
and qpa,o = ipi,a are imaginary. The component qt.hp is 
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always real as stated above. With these specifications we 
can rewrite Eqs. 


Wo 


= Nsiqt 


Sgji cos cosh (piz) + - qf) 

(10a) 


X cosh ( ) cos{zqt) 


c'qii 


= N. 


.911 


2qtpi cos 


X cosh ( ^ ) sin(zg() 


sinh(piz) + (g2 _ 

^ (10b) 


I ,a,i',x — ^aPt 


29 y sinh sinh(p,z) - (p^ + 

(10c) 


X sinh ( ) sinh(2:pt) 


= -iN, 


9|| 


2ptpi sinh ) cosh(piz) 

PiL 


-{Pt + 9 ||) sinh — cosh( 2 :pt) 


(lOd) 


In the same limit, the normalization constants are 

= V, for m = 0 and = 14/2, for m > 0, 


N, = Ai^Aqtq(\COs\qtL/2) 

,L 




-{pi-ir, 


L 


+ (9t - 9y) cosh (piL/2.) 


2 2^sin(9ti) 


+49y9t(9? - qf) cosh^(piL/2) sin(gtd) ^ (11a) 


N„ = ^ \ ^Pt^ii sinh (jptL/2) 
L 


, 2 , 2Nsinh(pzL) 


+(p? - «i)f 


+ (Pt + 9 y) sinh (jpiL/2) 


, 2 , 2Nsinh(ptL) 2 

-49^Pf{P( +( 7 ^)smh^(p,L/ 2 )sinh(p,L) lllb) 


C. The interaction 


Electron-phonon interaction in the volume of the film 
is described by the deformation potential HamiltoniarP^, 


Hdef = -EF d\¥{v)M>{v)V-uiv), (12) 

O Jv^i=Axd 


where Ep is the Fermi energy. From Eq. (12) we notice 


that all the transversal components of the h, s, and a 


modes are not contributing to the electron-phonon inter¬ 
action. Plugging Eqs. Q and ^ into (12) we get 


Hdef - 

k|| .qn ,i,n,n' 


n' ,n t 

55,q||Ck||-|-q||.n'Ck||,naq||{ 


/ n\n\ 


-qyi 


I Ckii ,n^, 


t 

qiiC 


(13) 


where n and v! in the notation o? refer to the electrons 
states (see Eqs. [^and[^. The coupling constants are 

, 9 I h 


%q = -^EpNo,, 


2pWq||,^ JL/2-d 

9Wq||.^,z(z) 


*C1|1 ■ Wq||,5(2;) -f 


dz 


K'iz)4’niz) 

dz. (14) 


Because the h modes do not contribute to the interac¬ 
tion (the divergence of any transversal displacement field 
is zero) and we are interested in the low temperature 
limit, in the followi ng w e shall take into account in the 
expression of 9^q" (14) only the modes with ^ = (s, 0) 
and (a, 0). 

The energy transferred from electrons to phonons in a 
unit of time (heat flux) is 


P 


2 ^ SWqii ,5 

k||;k',q|| 


-pem 
^ ^,n,n 


'(k|i;k||,qii) 


■pab 
^ ^,n,n 


'(k|i;kl|,q„) , 


(15) 


where the factor 2 comes from the electrons spin degen¬ 
eracy. 

As we wrote in the introduction, due to the weak cou¬ 
pling between electrons and phonons we may assume that 
the electron system has an equilibrium Fermi distribution 
corresponding to a temperature Te and the phonons have 
a Bose distribution of effective temperature Tph- We shall 
use the notations /3e = l/^ksTo) and jdph = {ksTph), 
and (C may be omitted). With these no¬ 

tations and applying the Fermi’s golden rule, we obtain 
the emission and absorption rates T®™ and 


r5“x(k|i; 

k||- 

-qii,qii) = 

27r 

h 

1 n' ,n 

1 [n-iPph^qii 

)+l] 

^ / (/3e£fc|| ,r, 

c)[l 

- /(/SeCkii- 

-qiib 

.)] 



X(5(ek||-q|| 

,n 

ek||,n' + Cgy 

), 



(16a) 

r|!;.n'(kii 

-qii 

;k||,qi|) = 

27r 

h 

1 n' ,n 
^C.qii 

1 ^(/?p/t£ij||) 


X 

1 


n')]/(Xek||- 

-qy 

.n) 



X<^(ek||-q|| 

,n 

ek||,n' + e(j|| 

)• 



(16b) 


where f{x) and n{x) are the Fermi and Bose distribu¬ 
tions, respectively. Using the identity f{x)[l — f{x—y)] = 
n{y)\f{x — y) — fix)] and taking into account the S func- 
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tions in Eq. we can write 

(k|l; kj| - qy, qy) - (ky - qy; ky, qy) 

27r, 


= ” /(/3eefe||x)] 

X [?^(/3ee9ii) -)] (17) 

and the power flux becomes 

P = P^°^T,)-P<-^\T,,Tph). (18) 

D. Long wavelength approximation 

We find the relevant low temperature asymptotic ex¬ 
pressions by expanding all the quantities in a Taylor se¬ 
ries in gy. From Eqs. (|^ and Q we derive expressions 
for qts, Pis, Pta, and pia, which we then use to calculate 
all the other quantities. For the symmetric modes it is 
sufficient to express qt^s and pi^s to the lowest order in gy, 
whereas for the antisymmetric modes we have to express 
Pt^a and pi^a up to the third order. In this way we obtain 


Qt,s 

Pl,s 


N, 


Pi,a 


OJ 

-2 


gy\/3 — 4J, (19a) 

gy(l-2J) (19b) 

2gyCtVl - J = 2gyCi-v/j(l - J), (19c) 

16g[|(3-4J)(l-J)2, (19d) 


9ii |i “ ~ + 


2 , 2J(1 - J) 


27 - 20J 


- ^(1 - J) 




(20a) 


Pt,a 


9|n 1 “ w(l “ '7)9|| + 


2(1-J) 


27 - 20J) 


- (1 - J) 


L 


CtL\ 


/l - d 9 

r ,1 

V 3 



N- 


9|| 


12_^6 


36 


9y \ (20b) 

(20c) 
(20d) 


Plugging Eqs. (10), (11), (19), and (20) into (14) we 


may write in the first stage 

n',ra|2 _ 


= Wsyhq\\{q'ii-pf)y^^^cos^{Lqt/2) 
xS{n',n,s,pi), (21a) 

KiP = \Nay[ptqiliq^-pf)f^smh\Lp,/2) 

xS{n' ,n,a,pi), (21b) 


where the overlap integral S{n'^ n, a^pi) is in cases s and 
a 

2 


S{n',n,s,pi) = 


S{n',n,a,pi) = 


hji 

J 

L/2-d 

L/2 

/ 

L/2-d 


(j)l, (z)(j)n{z) cosh(zpi) dz 


(22a) 


(j)*^,{z)(j)n{z) sinh(zpi) dz 


(22b) 

Using Eqs. ([^ and (19) and expressing all the qu antit ies 
In terms of the frequency uj we obtain from Eq. (22a) in 
the first relevant order 


/S'(ni, 712, S, w) « 1, if Til — 712 = 0 


(23a) 


^d^ (1 - 2J)^ 


1 


n 2 


(771-772)2 (771-1-772)2 


IGTT^cf J2(l _ J)2 
if 77i — 772 = 2777 
_ UJ^{L-d/ 2 fd? ( 1 - 2 J )4 

J '^{1 — jy \_{ni — n2Y ( 771 - 1 - 772 ) 2 ] ’ 

if 77i — 772 = 2777 -|- 1, (23c) 


1 


1 


(23b) 

n 2 


whereas from (22b) we get 


o, , _ W 1 (7 - <i)" 

if 77l — 772 = 0 

^ 9v^ d'^{L-df 3 
^ L^C^l[J{l - J)] 3 / 2 "' 

if 77i — 772 = 2777 

4^3 


w 


1 


1 


(24a) 
1 2 


-LO 


(771 - 772)2 (77i-h 772)2j 

(24b) 

1 1 " 


a/J(1 - J)Lci [(’21 - 772)2 (771 -I- 772)2j 

if 77i — 772 = 2777 -|- 1 (24c) 

Using these approximations we can now calculate the low 
temperature heat flux. 

III. ELECTRON-PHONON HEAT FLUX 

We separate and P^^'> into the symmetric and the 
antisymmetric parts, according to which type of phonons 
contribute to the process: P^°^ = P^''’ + pP and P^^^ = 
Pi^^ + P^^y We work in the low temperature limit, using 
the approximations of Section [IID[ 

A. Calculation of 

1. Contribution of the s modes 


Using (19) and (21) and changing from summations 


over ky and qy to integrals, for the symmetric part we 
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have 


integrals taken separately from Eq. (25), so 


a2 _ _ roo poo p2TT 

^ ~ ^(27r)3 EEX dkw fc|| J dqi\qii j d(j) 

xS'(n',n,s,q;)fia;5||n(/3ee,||){/[^e(efe|| - £ 9 , 1 )] 

~/(i®eefc|| )}ij(£k||-q|| ~ ^k|| + £g|| ) 

2A hEl J v-v- r „ , /■“ , 

= <‘kkl ‘i«ii9ii 


p27T 


d(l)Uj'^S{n',n,s,qi)n{l3eeph){f[Pe{(^k - ^ph)] 
-/(,5eefc)}<^(ek-qii - Ek + ^s) (25) 

The energy conservation gives the angles (j) (which may 
take two values): 

^ ^ + k?{kl + ql-2k\\q\\ cos 4> + kl 2 ) 

2me 2TOe 

—hujg (26) 


pOO pOO 

5 ^°^ = 2 ci^ 2 meJ{l - J) J dcph J 

X , {/[/3e(£e - Cp/i)] - fUde^e)} 

€e - e)^} {n\ eph) 


_ pOO 

= 2ci^2m,J(\ - J){kBTey^^ V / dxp,, 

/ -/O 


^ph 




dx 


1 (y ^ 71 ^ (^^i^phi^e)) ^ph _j_ 


^x-{y-x''^} {n' ,Xph,T^)) _|_ ^ 


where Xp/i = /3eep/i, iCe = /^eCe, a:„/ = 

xl^}{n,Xpk,T^} = PeS^hnAph), V = Pef^F, X -- 

xl^}{n,Xph), and 


(30) 


jde^n' 5 

: Xp - 


We change the variables in Eq. (25) into electron and 


phonon energies and denote them by Cph and Ce, respec¬ 
tively. Then Ce takes values from e„/ to 00 . Nevertheless, 
the energy conservation imposes extra constraints on the 
electron energy, such that the lower limit increases from 
£„/ to 


y - xl^}{n', Xph, Te) = y- Xn' - 


2mecfJ(l — J) 


ksTexl^ 


Xph + 


ksTeXh^ 
8mecfJ{l — J) 


(31) 


(«2, eph) = Eni + (^2, eph) (27) 


where 


Sei^}in2,eph) = 


2me J(1 — J)cf 


[iph 


^ph 


^ph 


mecf8J{l — J) 


(28) 


(s) 

in general we shall use the shorthand notations eli/ and 
(5ei*\ when this does not lead to confusions. 


Taking the three cases (231 for v! and n we write 
as the sum of three terms: 


From Eqs. (30) and (31) we calculate pfJiTe), which 
is plotted in Fig. ^as a function of Tg and d. We observe 
very sharp crests forming in narrow d intervals, which 
correspond to Ep « e„y, as explained in Appendix [X) In 
in the regions between the crests 


P) 


(0) 


Is 


{keTef 


E 


J 


TT^ me^^E% 


135^2 h^Lpcf Y (1 - J)^ 


^ ^/Ep — 


(32) 


Pi°) = 


rUeEp 1 

86Tr'^Lpcfh^ (1 — J)^ 

(0) I z:){0) I pi^) 




+ sS’+S®') 


_ p\^J I p 

— ^Is ' ^2s 


3s 


(29) 


The terms S^s ^ and are the summations and 


whereas on the crests the heat flux may increase by more 
than one order of magnitude and is not proportional, in 
general, to a single power of Tg. 

The contribution coming from the electrons which are 
scattered between sub-bands of indexes that differ by an 
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FIG. 2: PisIVei (Eqs. 29 and 301 vs Te and d for L = 100 nm. 
The lower plot is a zoom in the upper figure. forms 

narrow crests approximatively in the Tf, direction. 


where {n,Xph,Te) and 

(s) Q 

y - ^n' = y-Xn' - P, 


(s) _ n 2mecfj{l - J) 


^ph 




(34) 


For we have 


n —n—2k-\-l «oo 


^'is = X! / / dp€lf^S{n,n,s,qi) 

n',n ® €i 0 

Xn{Peeph){f[Peiee - eph)] - fiPe^e)} 

^ ^(^k|| —qy ,n Ck|j ,n' T ^pfC) 

{L — d/2)'^cP‘yj2me (1 — 2J)^ rp', 15/2 

27r4S4cf [J(l-J)]3/2^ ^ 


n^—n2=2fc+l n 


1 


E E 

n' 712 

r°° dxph^ r°° dxe 

Jq e^pf' -1 Jo ypx 

1 


1 


n 2 


{n'-n 2 Y (n'+ n 2 Y 

1 




^x-{y-x^^/)-Xph _j_ 2 ^ 


(35) 


where again 2 ;^ = xJJ{n,Xph,Te). sJJ and sJJ are 
plotted in Fig. (top and bottom, respectively). Com¬ 
paring these plots with Fig. [^we see that 5*2°^ and sJJ 
practically do not contribute to . 


even integer [n' — n = 2k) is (see Eq. 23b) 


Is' —n—2k 


-(0) 

^2s 


= E E 


de. 


ph 


dee 


dp e'. 


ph 


2 . Contribution of the a modes 
For the antisymmetric modes we have 

a2 _ _ poo poo p2'K 

" 2 (^EE "^^11 ^ 11 dp 


xS[n\n,s,qi)n{Peeph){f[Pe{ee - eph)] - /(/JeCe)} 

^ ^(Ck|| —qy ,n Ck|| ,n' F ^ph) 

_ 2J)4 n'-n=2k 


x|7V„| 


<7*911 


(<7y-iP?) 


2 16 ttE^ , 2 
sm 


9puja 


Lqt 

2 


87r4fi4c3 [J(l _ J)]3/2 


iksTr=/^Yl H 


1 2 


{n' — n)‘^ {n' + n)^ 


^Ih dxph 


dx 


1 


X S{n' , 122 , a, qi)hlOan{PeegY{f[l3ehk^^,n' “ ^ph)] 

,n') J’^(^k|| —qy ,n F 7 *LJq) 

2^4 irieEp J j 

p2-K 

X / (i(()S'(n',n,a,gi)n(/3eep/i){/[/?e(efc - Ep/i)] 


\/x y gX-iy-x'flY-Xph _|_ 2 ^ 

1 1 


-/(/ 3 eefc)}d(eky -qy ~ Ck F huJa) 


(36) 


gX-{y-x''Y) _|_ 2 


Plugging in the energy conservation from the <5 function, 
the power flux emitted from electrons to the antisymmet- 
(33) ric phonons can be divided in three parts, in the same 









































where the variables for the antisymmetric modes are 



= en' + 


-\-f-ph 


UleLci / J(1 — J) 1 


^ph 




2meLci y J(1 — J) 


= + Sel^,\n,eph), 


x^^)(n,Xph) = Pef^n'’{n,eph), and 


= R .(“) 


(n, Xph) = PeSe'^,’ (n, Cph) 


(a)/ 


rUeLci / J(1 - J) 

2hxph V 3 


Xn Xn' + Xpji 


h 


2meLci y J{1 — J) 


(39) 


We note that in the case of scattering on the antisym¬ 
metric modes 5x^^) {n^Xph) depends on Te only through 
Xnt Xji'^ and Xph. 

The function Pia/Vei is plotted in Fig. We notice 
similar crests as in Fig. which correspond to Ep « Cnp ■ 
The nature of these crests is discussed in Appendix 
In the regions between the crests the thermal power is 
described by 


FIG. 3: P 23 (top) and P 33 (bottom), for comparison with 
S[^J (Fig. and s[°J (Fig. [^. Both plots are calculated for 
L = 100 nm. 


way as we did before, 


p(0) 

A 


31/45 J3/4 

87l7r372^Uy (1^7)574 


n<nnp 

x(fcBTe)"/' 5 ] 

n 



(40) 


p(0) 

a 


2 A rUeEl J ( 0 ) ( 0 ) (o)^ 

3^2L3 (1 _ J) + *2a + -^3a 1 

Pia+P2a+P3a- (37) 


whereas on the crests the heat power may increase by 
more than one order of magnitude and is not propor¬ 
tional, in general, to a single power of T^. 


Using Eqs. (24b I and (24c) we write 


For the intra-band transitions we have 


(0) _ 1 hrn. 

^ T n 


1/4 


{L-d) 


1“ 2\ 2h^ [J{1-J)_ 

dxphx^j^ dx J 




gXph _ X Jp y/x 1 ^x-{y-x^y {n' ,Xph))-Xph _|_ ^ 



/■“ dxph Xpf 

dx J 

f 1 

/o e^p'* - 1 j 

0 Vx 1 

1 ('^’^pA)-^ph _j_ 


1 


^x-{y-x^^} (n' ,Xph)) _|_ 2 ^ 


(38) 


^x-(y-x*-^,\n,Xph)) _|_ 2 


(41) 
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FIG. 4: Pxajy&i (Eqs. |37| and [38| ) vs Te and d for L — 100 nm 
for L = 100 nm. The lower plot is a zoom in the upper figure. 
Pj® forms narrow crests approximatively in the Te direction. 


and 



1 




+ 1 


(42) 


do) 


The functions and are plotted in Fig. i (top 
and bottom, respectively) and we see that they do not 
practically contribute to the heat power. 

Adding all the contributions (29) and ( [37| ) we obtain 
the total heat power emitted from the electrons into 
the phonons system. is plotted in Fig. and the 

ratio is plotted in Fig. We see that the 

dominant contribution comes from the heat power to the 
antisymmetric modes. Theref ore i n the low temperature 
limit is described by Eq. (40) between the crests, so 



FIG. 5: P 2 s (top) and Pss (bottom) for comparison with 
(Fig. i and (Fig- |^- Both plots are calculated for 
L — 100 nm. 


it is proportional to Along the crests, both and 
may contribute substantially, together with Pa^\ 
and there is no simple temperature dependence of . 

For a rough comparison with experimental data pub¬ 
lished in Ref. mm we plot in Fig. I^the heat power for 
a sample with L = 45 nm (like sample Ml in Ref. [17)1. 
Although the experimental data is scattered, we observe 
a qualitative agreement. 


B. The heat flux from phonons to electrons 


The heat power P^^^ from phonon to electrons is cal¬ 
culated as P^^\ with the exception that in the integrals 
over eph one should use a Bose distribution corresponding 
to the temperature Tph instead of Tg. 






































1. The contribution of the s modes 


As in Sectionflll A 1[ we write 


= A 


rrieEp 


^Gn’^Lpcffi^ (1 — J)^ 




= P 


+ P^'J 


P 


( 1 ) 


3s ’ 




(1)^ 


(43) 


where 


_ _ poo poo 

5*^^^ = 2ci\/2meJ{l - J) / dcph / dee 

W do 

^Ih^iPph^ph) 


= {/[/3e(ee - Cp/i)] - /(/3eee)} 


\Jee- e^hn'^f^ph) 

p OO 

= 2ci^2meJ{l- J){kBTey'^Y. I dxph^^£d_ 

m' ^ 


1 


dx I 

y/x 1 ^x-(y-xy (n',Xph.,T^))-Xph _j_ ^ 

1 1 


^x-(y-xy {n',Xph,Tc)) _|_ 


(44) 


where t = Tg/Tph, and all the other quantities are the 
same as in Se etio n riII A II 

From Eqs. (44) we calculate Py\Te), which is plotted 
in Fig. as a function of Tg and d. We observe very 
sharp crests forming in narrow d intervals, which corre¬ 
spond to Ep Ki Cnp, as explained in Appendix [C| In the 
regions between the crests 




T [K] 


FIG. 7: The ratio for ^ = 100 nm. We see that the 

dominant contribution comes from the heat flux from elec¬ 
trons to the antisymmetric phonons. 





FIG. 8: (Color online) p(°) as a function of Te for different 
values of d (shown in the legend) and L — 100 nm. The 
higher values of are from the crests regions (see Fig. [^. 
We notice that the heat flux cannot be written as xTf. 


whereas on the crests the heat power may increase by 
more than one order of magnitude and is not described 
by a simple power law. 


FIG. 6: p (°)/Vei for L = 100 nm. 


The contribution coming from the electrons which are 
scattered between sub-bands of indexes that differ by an 
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even integer {n' — n = 2k) is (see Eq. 23b) 


2 ' —n—2k 


= E E 


n' n 


I dcph I 
lo Jei 


dCe 


d4>^lh 


xS{n',n,s,qi)n[l5ph^ph){f[Pe{^e - e-ph)] - f[Pe(^e)) 
x<5(ek|| — q||,n ^k|| ,n^ H“ ^ph) 

>■4 n' —n—2k 


V^d^ (1 - 2J)^ 

STT^h^cf [J(l- J)]3/2 




1 


1 


n 2 


/O 


{n' — n)^ (n' + n)^ 

r°° dx \ 1 

lo 1 _|_ 

1 1 


n' n 

xlh dxph 

^tXph _ ^ 





whereas the rest of the terms in the summation (43) are 

n' — n—2k-\-l „qq „qq „2-k 

= E / d€ph dce d(j)€l^S{n',n,s,qi) 

n'.n do Je, Jo 

Xn(/3p?,ep/i){/[/3e{ee - Cph)] - f{Pe^e)) 

^ —qy ,n ^k||,n' “t“ ^ph) 

{L-d/2)^d^y/2^ (1-2J)4 


2Tr*h'^cf [J(l - J)]3/2 

n' — n 2 — 2 k-\-l 

xE E 


1 




(n' — n-iY (n' + n 2 )‘^ 


^Ph /■“ dXe 


/g — 1 7o 


x-{y-x^^})-Xph _j_ 2^ 


pX-(y-x^‘}) 2 ^ 



FIG. 10: Pi'.^VVei (Eqs. _ _ 

plot is a zoom in the upper figure. 


29 


and 


301 vs Te and d. The lower 


Like in the case of the terms ^ 2 ^^ do not 

contribute significantly to pI^\ 


3s 



14 0.01 j [K] 

e 

FIG. 9: P^°\T,,d)/V,i for L = 45 nm. 


2. The contribution of the a modes 


For the antisymmetric modes we have 


pI^'> = 


2A rrieEp J 


Stt^L^ pc'jfo (1 _ j) 


is[l^ + + sg^) 


(1) , d(1) 


— 4^1a T X2a 


P- 


( 1 ) 


3a ’ 


(46) 


where 


(1) ^ 1 ^ 
2\ 2h? 


J(1 - J) 


1/4 


(L-d) 




dxphX^J^ r°° dx 


1 


^x-(y-x)f,\n' ,Xph)) _|_ 2 ^ 


0 Vx 1 gX-(y-x''fff (n',Xph))-Xph _|_ ^ 

(47) 
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cn 



On 


FIG. 11: Pia'/Vei VS Tf, and d. The lower plot is a zoom in 
the upper figure. 



FIG. 12: Total power flux per unit volume jV^i from the 
phonon system into the electron system for L = 100 nm. 


The other summations of Eq. (461 are 


c(l) _ 

^1a — 


9 

j me d‘^{L — dY 

3 

27r4\ 

2fi7 (Pci)5/2 J(1 - J) 



1/4 


n' — n—2k 

X E 


n' ,n 


1 1 
{n' — nY {n' + n)^ 


{ksT^f 


dxphX^Jh dx 



g_x-{y-x'i) _|_ 2 


and 


The variables are the same as before. 


The function PialVei is plotted in Fig. 11 and we no¬ 
tice similar crests as in Fig. [T0| T he nature of these crests 
are studied in the Appendix |D| In the regions between 
the crests 


P. 


( 1 ) 


la 


{kBTphf'^ 


31/45 /7 

8 ^ 2713/2 ^ [2 


ml^^Ep{L - dY 
h^/^pc/^LV^ 


J3/4 

X(i_j)5/4 ^ 


1 

Y Ep — < 


(48) 


whereas on the crests the heat power may increase by 
more than one order of magnitude and is not described 
by a simple power law. 


S. 


( 1 ) 

3a 



/■“ dxph xH^ 

dx 

1 

Jq — 1 

Jo Vx 

^x-{y-x[)-Xph 


1 


^x-(y-x{) _|_ 


+ 1 
(50) 


but they do not contribute significantly to the total heat 
power. 

Adding all the contributions to P'd-'> we obtain the total 
heat power from t he phonons to electrons system. P^^'^ 


is plotted in Fig. 
in Fig. 


13 


12 


and the ratio is plotted 

that has the 


We observe in Fig. 


13 


dominant contribution and th erefore between the crests 
is described by Eq. (48), with a temperature de¬ 
pendence oc On the crests there is a combined 
contribution of both, Pi^^ and P^\ 
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FIG. 13: The ratio /Ps^^ for L — 100 nm. We see that 
the dominant contribution comes from the heat flux from elec¬ 
trons to the antisymmetric phonons. 


where s is the smaller of the dimensions of the electron 
and phonon gases.^^ On the other hand, in an experi¬ 
mental setup various factors, such surface contamination, 
thickness variations, or other imperfections might make 
it difficult to observe the strong oscillations of P with the 
thickness of the film and eventually an averaging proce¬ 
dure might be more appropriate. 
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IV. CONCLUSIONS 


Appendix A: The contribution 


We calculated the heat flux P between electrons and 
phonons in a system which consists of a metallic film of 
thickness d of the order of 10 nm deposited on an in¬ 
sulating SiNj, membrane of thickness L of the order of 
100 nm (see Fig. [^. We described the electrons as a gas 
of free fermions confined in the metallic film. The prin¬ 
cipal characteristic of the electron gas is that due to the 
quantization of the wavevector component perpendicu¬ 
lar to the film surfaces the electron gas forms quasi 2D 
sub-bands , which ar e specific to a quantum well (QW) 

formalism .CSElEIEn] accompanying paper we de¬ 

scribe the electrons as a 3D Fermi gas. We choose a 
range of parameters L, Tg (electrons temperature), and 
Tph (lattice temperature) such that the phonon gas has 
a quasi two-dimensional distribution (for L = 100 nm we 
took Te, Tph < 0.2 K). 

The heat flux assumes the non-symmetric form P = 
p(o)(Te) — {Te,Tph), where is the heat flux 

from electrons to phonons and P^^'^ is the heat flux from 
phonons to electrons. Due to the sub-bands formation, P 
exhibits very sharp crests approximately along the tem¬ 
perature axis in a |P| vs (d, T^) or |P| vs (d, Tph) plot, sep¬ 
arated by much more flat valleys. In the valley regions, 
in the low temperature limit P oc , whereas 

on the crests the heat power increases by more than one 
order of magnitude and P does not obey a simple power 
law behavior. Such a sharp oscillatory behavior is a char¬ 
acteristic of the Q W descript ion of the electronic states 
in the metallic filmalSElBIIISI ^nd could be useful for ap¬ 
plications, like thickness variation detection or surface 
contamination. Moreover, by varying the thickness of 
the membrane one might vary dramatically the electron- 
phonon coupling and, through this, the thermalisation or 
the noise level in the system. 

The power-law behavior P oc — does not agree 
with the generally expected formula P oc 


From Eq. (29) we have 


p(0) _ A _ 


mgE 


e^F 


4°) (Eq. 

29). 

c.(0) 

(Al) 


where 

5i°^(y,Te) = 2Qy/2med(l- J)(fcBTe)"/" 


(A2) 




r°° dxph xlf^ dx 


1 


— 1 do \ \ 

^ ' (A3) 




and Sxn'^ = Sxn^ (n,Xph,Tg). We take separately the 
integral over x and generalize it to 


h{r,y) = [ dxx'^ ^ (- 
Jo \( 


1 


x-(y-5xi!’'')-Xph _|_ 1 


1 


v-{y-Sxi^^) _|_ I 


(A4) 


If y — 6xn' ^ 1 in Eq. (A4 ) we can use the approximation 

r — 1 




(AS) 


Taking r = 1/2 and plugging the highest order approxi¬ 
mation back into (A3), we get 


1 

I„(o}{y,Te) ^ —p I dx. 

*1” Vy JO 


‘'ph 


TT^ 1 

= - —. (A6) 


(e^p'* — 1) 15 ^ 
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FIG. 14: (Color online) (a) Umax vs Te, (b) hs{y max ) vs. Te 
and (c) —J/max(re)] for Te = 0.1 K (blue curve) and Te = 

0.2 K (red curve); the two curves in (c) are indistinguishable. 


The approximation (A61 cannot hold it y — Sxn' is of 
the order 1. In such case we have to evaluate Eq. (A3) 
explicitly and we notice that Ig(o){y,Te) forms a maxi¬ 
mum at y = ?/max('7e) (see Fig. [13. The shape of /^(o) as 
a function of y — j/max varies little with (Fig. 14 c) and 
the maximum value /g(o)(ymax) ~ 4.5 slightly decreases 

with Te (Fig. [I^b). 

In the limit y <C —1, /„(o)(y,Te) decays exponentially 

^Is 

with y and therefore we neglect all terms with n > np + t 
from the summation (A2). The summation extends up 
to np or np + 1 and we may have three situations: 


S 


(0) 

Is 


{ 4 riF-l ^ 

}-{kBTef Y. rr - 

15 VEp-en 

-h(fcBTe)^/^/g(0)(y - x„y)|, (A 7 a) 


when y - Xnp - 


> 0 and the approximation (A6) 


does not hold for I^(o) (y — Xnjr,Te), 

'^Is 


0^4 1 

S[°J « ClV^meJ{l-J){kBTerY rr - - 

J-b V-Ff - e-n 

(A7b) 

(s) 

when y — Xnp — Sxnp 1 and y — <C —1, and 

finally, 


^(0) 

’is 


2ci^2meJ{l-J) -p{kBTef Y m - 

I 15 ^/Ep - e„ 

T{kBTe)'^I^I^(Yy ^np + l )}, (A7c) 


when y - Xn^+i < 0 and /^{o) (y - Xn^+i^Te) is not neg¬ 
ligible as compared to the whole summation in (A21. 

We say that in cases (A7a) and ( ATcb we are in the 
crests regions of whereas in case (A7b) Pis cx Tg and 


the heat power may be more than one order of magnitude 
smaller than in cases (A7a) and (A7c). 


Appendix B: The contribution Si°^ 


From Eqs. (|37|) and (|38|) we write 

J 


p(o) _ 

-^la 


2A 


meEj, 


(0) 


Stt^L^ pcjJA (1 — J) 


(Bl) 



2h^[j{l-J)\ 


1/4 


(L-d)^ 


{keTef 






(B2) 


where 

isf) {y) 


r°° dxph xl^^ 

J g^^ph _ X 

1 


r— 

'o 


^x-(y-5x\^^)-Xph _j_ 2 

(B3) 


and x''n’ = xit\n,Xph) (Eq. 

We observe that the integral over x is similar to Eq. 


39). 


A4 


In this case Sx''n'^ depends on the variables n and Xphj 
but not explicitly on Tg, like Sxn \n,Xph,Te). Therefore 
1^(0} (y) has a shape which is independent of Tg (see Fig. 

[isl with a maximum at ymax ~ 5.6. 

If y 3> ymax we can use the approximation (A5) to 
write 


IgW (y) 


\/y^ 


dxph xH^ 

gXph _ 1 



(B4) 


where C(3-5) ~ 1.13 is the Riemann’s zeta function. 

If y <C —1, 7„(o)(y) decays exponentially with y and 

*10 

therefore the terms /„(o)(y — a:l“^) with n > np + 1 will 
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when y — Xnp+i < 0 and Ig(o) {y — Xnp+i,Te) is not neg¬ 


ligible as compared to the whole summation in (B2). 


In the cases ( B5a ) and (B5c), d cor responds to crests 
regions of Fig. whereas in case (B5b), d corresponds to 
the valley region and the heat power may be more than 
one order of magnitude smaller than in the crests. 


Appendix C: The contribution 


( 1 ) 


From Eqs. (43) and (44) we have 


^ (Cl) 


FIG. 15: I„(o) vs y has a maximum at y Ri 5.6. /„(o) (y') 

*l<i ^la 

decays exponentially for y —1 and is proportional to 
for . 


mn^Lpcfh^ (1 - J)2 1* ’ 

20 where 

= 2ci^/2m,J{l-J){kBT,y^^ 

X ^ Ig(iJ {y- Xn,Te, Te/TphlP2) 


be omitted from the summation (B3). The summation 
extends up t o np or and we 

like in Eqs. (A7): 


rave again three cases, 


S- 


(0) 


1 [rrie 


1/4 


[L-df 


2V 2h^ [j(l - J)j 

X \ ^gW {y ^np) 


(B5a) 




\/ Ep — I 


when y — Xnp — Sxl^l > 0 and the approximation (B4) 
does not hold for J„(o) (y — Tg), 

*la 


^(0) 

’la 


16 


-cU 


n<n. 




'We 

3 

2h^ 

[j(l-J) 

r.p 

~y 

I 


3 1 (L - df 


y/Ep — I 


(B5b) 


when y — Xnp — dx^J ^ 1 and y — Xnp+i <C —1, and 
finally, 



(B5c) 




n<njr 


y/ Ep — 


r ! rp r r 

Jo e^-.'^-lio 


0 V® 
1 


QX-{y-Sxi^^)-Xph _|_ \ gX-(y-Sxi‘'^) _|_ ^ 


(C3) 


with similar notations as in Appendix [Aland t = Te/Tph- 
From Eq. (C3) and using the results ot Appendix we 


(s) 

obtain in the limit y — 6xn 1, 


Ig(r,{y,Te,t) « 


y/y. 


dx. 


‘'ph 


TT^ 1 


P^ l^gtxph — 1 ) 15 

(C4) 

If y in Eq. (C3) is of the order 1 (positive or negative), 


then we evaluate Eq. (C3) explicitly. Taking Tph = 0.2, 
/g,(i) (y, Te, Te/Tp/i) forms a maximum at y — ymax(Te) 
(see Fig. [l^ . 

In the limit y < -1, I^m{y,Te,Tf,/Tph) decays expo¬ 
nentially with y and therefore we neglect all terms with 
n > Up + 1 from the summation ( |A2[ ). 

The summation extends up to or np + 1 and we 
may have three situations: 


di) 


^^{kBTpuf 

n^ —1 

E 

n—1 

Xnp) 1* ; 



1 


— Cn 

(C5a) 


when y — Xnp — dxnp ^ 0 and the approximation (A6) 


does not hold for {y — Xn,=-.Te 
Fi) 


« —Ciyj2meJ{l - J){kBTpnfY.^ 

n—1 ^ 
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(C5b) 
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FIG. 16: (Color online) (a) ymax vs Te and (b) 7^(i) (t/max) 
vs. Te for Tph = 0.2 K. 


(s') 

when y — Xnp — Sxnp ^ 1 and y — Xnp+i ^ —1, and 
finally, 


S 


( 1 ) 

Is 


{ 4 np 

— {ksTph) , — 

„=1 V^F - en 

+ {kBTef^'^Ig{i){y - Xnp + l)^ , (C5c) 


Appendix D: The contribution S 


( 1 ) 


From Eqs. (46) and (47) we write 


p(i) _ 

J 1 ^ — 


c(i) _ 

“ 2 V TqX 


2A irieEp J 

pcfh^ (1 — J) 


c.(l) 


(Dl) 


1 [me 


1/4 


(L-d)^ 


{keTef 

y.Y,Ig^.){y-x^:\TjTpH). (D2) 


where 


r 1 r 

dg(^){yji) I tx h 1 I r~ 

* 1 “ Jo “ 1 Jo 


1 


gX-(y-Sxt^) _|_ I 


'Jx \ QX-(y-5x^^^)-Xph _|_ \ 

(D3) 


and keeping the rest of the notations like in Appendices 
EandO 

has a maximum at y = ymax, which is plot- 


•SL 

ted in Fig. j^(a) together with /^p) (j/mL, t) (Fig.|l7|b). 
If y 3> 2/max we can use the approximation ( |A5[ ) to write 

1 dXphxH^ _i5^rr\ 






(D4) 

If 2 / 'C —1, Ia(i)(y,t) decays exponentially and therefore 

“^la 

the terms /„(i) (y — x^^) with n > np + 1 will be omitted 

“^la _ 

from the summation (D3). 

The summation extends up to np or np + 1 and we 
have again three cases, like in Eqs. (A7): 



n<nn 

+(tBT,ov>i^c(n ^ 

\ / ^ 


(D5a) 


y/Ep — 


when y — Xnp — 6xnp ^ 0 and the approximation (D4) 
does not hold for J„(i )(y — Xnp,T^), 

^Is 


when y - Xnp+i < 0 and IgO) {y - Xnp+i,Tf,) is not neg¬ 


ligible as compared to the whole summation in (A2). 


Equations (C5a) and (C5c) describe the crests regions 


of whereas in case (|C5b| Pi,, (x Pjf, and the heat 


wnereas m case ^ 

power may be more than one order of magnitude smaller 
than in cases (C5a) and (C5c). 
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finally, 





(D5c) 




y/ Ep — 


when y - Xnp+i < 0 and /^(o) (y - Xnp+i,Te) is not neg¬ 
ligible as compared to the whole summation in (B2). 


FIG. 17: (Color online) (a) ymax vs Te and (b) /^(i) (t/n 
vs. Te for Tph = 0.2 K. 


In the cases ( D5a ) and (D5c), d cor responds to crests 
regions of Fig. whereas in case (D5b I, d corresponds to 
the valley region and the heat power may be more than 
one order of magnitude smaller than in the crests. 
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